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We  investigate  CO2  bubble  dynamics  on  the  anode  side  of  a  direct  methanol  fuel  cell  (DMFC).  In  con¬ 
trast  to  previous  studies,  we  analyse  the  effect  of  both  channel  wall  and  diffusion  layer  wettability  by 
observing  two-phase  flow  from  the  side  at  different  mean  velocities  of  the  fuel  supply.  Hydrophobic  and 
hydrophilic  flow  channel  surfaces  are  compared  experimentally.  The  hydrophilic  flow  channel  leads  to  a 
minimum  pressure  drop  along  the  channel.  Bubbles  show  virtually  no  pinning  and  consequently  travel  at 
approximately  the  mean  fuel  velocity  inside  the  channel.  In  contrast  to  this,  we  observe  bubble  pinning 
in  the  hydrophobic  flow  channels.  The  critical  fuel  velocities  necessary  for  detachment  of  the  bubbles 
mainly  depends  on  bubble  length.  We  identify  and  describe  a  new  bubble  bypass  configuration  where 
fuel  bypass  channels  are  solely  generated  in  a  favourable  position  underneath  a  blocking  bubble  along 
the  diffusion  layer.  This  enforces  fuel  to  bypass  the  C02  bubble  at  a  large  relative  velocity  close  to  the 
diffusion  layer,  thus  enhancing  mass  transfer.  Our  experimental  findings  are  in  excellent  agreement  with 
a  CFD/analytical  model.  This  model  allows  for  quantitative  prediction  of  average  bypass  flow  velocity. 

©  201 1  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Comparing  hydrogen-fed  polymer  electrolyte  membrane  fuel 
cells  (PEMFCs)  with  direct  methanol  fuel  cells  (DMFCs),  the  latter 
offer  some  unique  advantages,  that  have  recently  led  to  commer¬ 
cially  available  systems  focusing  on  off-grid  applications.  The  most 
notable  gain  is  the  simpler  fuel  storage  and  refilling,  as  the  liquid 
methanol  solution  employed  in  DMFCs  is  much  more  convenient 
to  handle  than  highly  volatile  hydrogen  [1  ]. 

Nevertheless,  there  are  also  disadvantages.  While  PEMFCs 
employ  a  pure  gas  with  ideally  no  waste  production  on  the  anode 
side,  the  methanol  used  in  DMFCs  must  be  diluted  to  a  concentra¬ 
tion  of  approximately  one  molar  methanol  in  water  to  mitigate  the 
losses  from  methanol  crossover  [2].  Additionally,  C02  is  produced 
as  a  waste  gas,  creating  two-phase  flow  including  blocking  bub¬ 
bles,  with  the  void  fraction  increasing  along  channel  length,  as  first 
observed  by  Argyropoulos  et  al.  [3 ].  This  leads  to  a  much  more  com¬ 
plicated  anodic  mass  transfer  system  and  to  higher  mass  transfer 
losses  as  current  density  increases.  As  a  result,  the  current  density 
of  a  DMFC  is  lower  than  that  of  a  conventional  PEMFC,  due  to  higher 
losses  on  the  anode  side  [4]. 
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Actively  transporting  liquid  instead  of  pure  gas  through  the  flow 
field  of  a  fuel  cell  causes  additional  system  losses  due  to  an  increase 
of  necessary  pumping  energy.  The  pressure  drop  can  be  predicted 
and  optimised  only  by  considering  two-phase  flow,  as  the  pressure 
drop  induced  by  blocking  bubbles  due  to  interfacial  forces  between 
liquid  and  gas  is  significant  and  cannot  be  neglected  [5].  Therefore 
understanding  bubble  behaviour  is  of  special  interest  for  design  and 
analysis  of  liquid  feed  direct  methanol  fuel  cells. 

A  common  method  for  analysing  two-phase  flow  in  fuel  cells  is 
visualisation  in  a  transparent  flow  field  which  allows  a  top  view  of 
parts  or  the  whole  flow  field  [3,6-8].  An  alternative  method  pre¬ 
sented  in  [9]  detects  void  fraction,  slug  velocity  and  liquid  slug 
length  at  the  exit  of  a  DMFC  anode  flow  field.  Similar  techniques 
are  also  employed  in  general  microfluidics  [10-15].  Fuerstman  et  al. 
[16]  derive  the  pressure  drop  in  a  rectangular  micro-channel  with 
flowing  liquid  and  blocking  bubbles  dependent  on  channel  aspect 
ratio  and  the  liquid  used.  Zhu  and  Petkovic-Duran  [17]  observe 
interface  velocity  when  employing  different  channel  materials  and 
surface  treatment  methods. 

Instead  of  regarding  pressure  losses  due  to  pinning  bubbles  in 
DMFCs  to  be  a  major  drawback,  some  groups  utilise  two-phase  flow 
to  circulate  fuel  in  a  completely  passive  manner.  They  employed 
buoyancy  [18]  or  tapered  channel  structures  [19,20]  to  supply  the 
fuel  to  DMFCs.  A  disadvantage  of  this  solution  is  the  lower  current 
density  due  to  decreased  fuel  velocity  in  the  flow  channels,  which 
increases  mass  transfer  losses. 


0378-7753 /$  -  see  front  matter  ©  2011  Elsevier  B.V.  All  rights  reserved, 
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Nomenclature 

A chan 

open  channel  cross-sectional  area  (m2) 

Acor 

bubble  bypass  cross-sectional  area  (m2) 

c 

correction  coefficient  for  flow  resistance  of  non¬ 
circular  channel  cross-sections 

Ca 

capillary  number 

d 

thickness  of  Toray  TP-060  diffusion  layer  (0.1 9  mm) 

F 

Faraday  constant  (F=  96,485  C  mol-1 ) 

hub 

bubble  height,  distance  between  diffusion  layer  and 
channel  floor  (m) 

hcor 

bubble  bypass  height  along  channel  wall  (m) 

Imax 

maximum  intrusion  of  diffusion  layer  material  into 
channel  (m) 

Imin 

minimum  intrusion  of  diffusion  layer  material  into 
channel  (m) 

I 

electric  current  (A) 

lv 

bypass  flow  rate  (m3  s-1 ) 

I< 

curvature  of  an  interface  between  two  different  flu¬ 
ids  (m-1) 

k 

in-plane  permeability  of  Toray  TP-060  diffusion 
layer  ( 1 .26  x  1 0-1 1  m2 ) 

l 

channel  length  touching  active  area  (m) 

Ibub 

bubble  length  from  tip  to  tip  in  flow  direction  (m) 

Icor 

bubble  bypass  length  (m) 

™CH3OH 

mass  flow  of  methanol  (kg  s-1 ) 

Mch3oh 

molar  mass  of  methanol  (kg  mol-1 ) 

P 

static  pressure  (Pa) 

Pc 

capillary  pressure  drop  over  an  interface  between 
two  different  fluids  (Pa) 

Pc,b 

capillary  pressure  drop  over  the  back  interface  of  a 
bubble  (Pa) 

Pcf 

capillary  pressure  drop  over  the  front  interface  of  a 
bubble  (Pa) 

A  p 

pressure  drop  (Pa) 

Apbub 

capillary  pressure  drop  over  a  blocking  bubble  (Pa) 

Apbyp 

capillary  pressure  drop  in  a  bypass  along  a  blocking 
bubble  (Pa) 

Qdl 

lateral  flow  rate  in  the  diffusion  layer  (m3  s-1 ) 

Qsp 

flow  rate  generated  by  the  syringe  pump  (m3  s-1 ) 

V 

velocity  vector 

vbub 

bubble  velocity  (m  s-1 ) 

Vchan 

average  channel  velocity  (m  s-1 ) 

Vcor 

average  bypass  fluid  velocity  (m  s-1 ) 

Wbub 

bubble  width,  =wchan  (m) 

wchan 

channel  width  (m) 

Z 

valency  number 

P 

dynamic  viscosity  (2-molar  methanol  solution: 
0.001212  Pas) 

0ad,DL 

advancing  contact  angle  of  2-molar  methanol  solu¬ 
tion  on  diffusion  layer  (°) 

®  ad, dyn 

advancing  dynamic  contact  angle  (°) 

@ad,stat 

advancing  static  contact  angle  (°) 

@ad,w 

advancing  contact  angle  of  2-molar  methanol  solu¬ 
tion  on  channel  wall  (°) 

@rec,DL 

receding  contact  angle  of  2-molar  methanol  solu¬ 
tion  on  diffusion  layer  (°) 

@rec,w 

receding  contact  angle  of  2-molar  methanol  solu¬ 
tion  on  channel  wall  (°) 

P 

density  (2-molar  methanol  solution:  988  kgm-3) 

O 

surface  tension  (2-molar  methanol  solution: 
0.059  Nm-1) 

subdomain  i 

8tti 

boundary  condition  i 

One  can  conclude  that  an  optimum  between  losses  due  to  pres¬ 
sure  drop  and  mass  transfer  losses  must  be  found  to  achieve  the 
maximum  efficiency  of  the  system.  With  this  goal  in  mind,  we 
study  bubble  dynamics  in  a  DMFC  and  employ  a  new  visualisation 
technique  by  observing  the  flow  field  from  the  side. 

In  contrast  to  previous  studies,  we  take  into  account  the  dif¬ 
ference  in  wettability  between  the  channel  wall  and  the  diffusion 
layer.  We  compare  a  native  cyclic  olefin  copolymer  (COC)  channel, 
which  is  hydrophobic,  with  a  COC  channel  coated  with  PDMAA-BP, 
which  is  hydrophilic.  We  vary  the  mean  fuel  velocity  and  observe 
pinning,  contact  angle  hysteresis,  bubble  velocity,  bubble  blocking 
and  bypassing  flow  to  gain  a  deeper  understanding  of  two-phase 
flow  and  its  consequences  on  pressure  drop  and  mass  transfer 
losses  in  DMFCs. 

The  experimental  work  is  supported  by  a  dynamic  two-phase 
model.  In  recent  years,  different  approaches  have  been  presented 
to  model  two-phase  flow  such  as  Lattice  Boltzman  [21],  dissipa¬ 
tive  particle  dynamics  (DPD)  [22],  molecular  dynamics  [23]  or  the 
standard  fluid  simulation  method  based  on  finite  volumes,  which 
is  also  referred  to  as  computational  fluid  dynamics  (CFD)  [24-26]. 
Especially  the  last  method  has  gained  considerable  popularity  in 
the  DMFC  community  due  to  its  capability  not  only  to  model  two- 
phase  flow,  but  also  to  include  interacting  mechanisms  such  as 
electrochemical  kinetics  or  heat  transfer  by  coupling  several  trans¬ 
port  equations  [27,28].  In  the  work  reported  here,  we  employ  a 
CFD/analytical  model  to  quantify  fuel  transport  in  the  bypasses 
along  a  blocking  bubble  in  the  channel  of  a  DMFC. 

2.  Experimental 

2.1.  DMFC  setup 

To  prevent  interference  by  neighbouring  channels,  a  single 
channel  flow  field  is  fabricated  of  COC.  The  side  wall  is  polished 
to  achieve  optical  transparency  and  the  resulting  side  view  of  the 
experiments  in  the  anode  flow  field  channel  is  documented  on  film 
using  a  magnifying  lens  and  a  p^Eye  camera.  This  flow  field  is  inte¬ 
grated  into  a  fuel  cell  as  illustrated  in  Fig.  1.  A  passive  cathode 
fabricated  of  SGL  porous  graphite  plate  is  used,  and  diffusion  lay¬ 
ers  (Toray  TP-060)  ensure  a  more  even  fuel  distribution  on  both 
the  anode  and  the  cathode  sides.  A  catalyst-coated  Nafion  115 
membrane  (active  area:  2.6  cm2;  anode:  3  mg  cm-2  PtRu;  cathode: 
1.3  mg  cm-2  Pt;  manufactured  by  BalticFuelCells)  provides  separa¬ 
tion  of  the  two  half-cells.  As  described  in  [20],  contacting  is  achieved 
by  a  100  p,m  thick  gold  foil  covering  land.  Two  endplates  made  of 
brass  on  the  cathode  side  and  a  double  layer  of  PMMA  and  alu¬ 
minium  on  the  anode  side  provide  contact  pressure.  The  anode 
endplate  is  designed  with  a  viewing  window  to  allow  an  addi¬ 
tional  top  view  of  the  flow  field.  Two  flow  channels  with  different 
wettability  properties  are  compared.  This  is  achieved  by  treating 
one  of  the  COC  flow  fields  with  PDMAA-BP,  which  creates  a  more 
hydrophilic  surface.  On  the  anode  side,  a  small  load  of  0.9  mA  mm-2 
induces  C02  production.  A  2-molar  methanol  solution  is  pumped 
back  and  forth  with  a  syringe  pump.  Over  time,  a  channel-blocking 
bubble  grows.  After  the  bubble  reaches  a  specific  size,  the  load  is 
turned  off  and  the  dynamic  behaviour  of  the  bubble  in  the  channel 
can  be  observed  at  different  fluid  velocities. 

2.2.  Microfluidics  of  two-phase  channel  flow 

The  average  velocity  of  methanol  solution  in  the  flow  channel  is 
derived  from  the  volumetric  flow  rate  of  the  syringe  pump  and  the 
channel  cross-sectional  area.  To  calculate  the  exact  channel  cross- 
sectional  area,  the  width  from  the  left  to  the  right  channel  wall 
and  the  height  from  the  floor  to  the  top  of  the  employed  channels 
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Fig.  1.  Schematic  diagram  of  experimental  setup:  (1)  reservoir  1,  (2)  syringe  pump,  (3)  DMFC,  (4)  transparent  (COC)  anode  (5)  side  view  of  flow  channel,  (6)  blocking  bubble 
(7)  membrane  and  diffusion  layers  (8)  cathode  (9)  reservoir  2. 
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Fig.  2.  Schematic  diagram  of  channel  cross-section  area  after  clamping,  with  diffu¬ 
sion  layer  intrusion  i.  The  channel  width  wchan  amounts  to  1.135  mm  (native  COC 
channel)  and  1.170  mm  (COC  channel  coated  with  PDMAA-BP),  the  channel  heights 
hchan  are  determined  to  be  0.880  mm  and  0.910  mm  respectively. 


are  measured  after  milling.  We  expect  the  pressure  of  clamping  to 
push  parts  of  the  diffusion  layer  into  the  channel  and  thus  reduce 
the  channel  cross-sectional  area  as  depicted  in  Fig.  2.  To  quan¬ 
tify  this  intrusion,  we  employ  a  Tesa-Visio  300  video  measuring 
machine  to  measure  the  intrusion  in  situ  through  the  top  viewing 
window  of  the  anode.  The  highest  diffusion  layer  intrusion  of  50  p,m 
is  measured  in  the  middle  of  the  channel  (Fig.  3).  At  the  edge  of  the 
channel,  where  the  diffusion  layer  reaches  land,  the  intrusion  is  still 
about  20  |xm  due  to  material  being  pressed  up  along  the  channel 
wall  as  can  be  observed  in  [29]. 


Fig.  3.  Diffusion  layer  material  intrusion  into  the  channel  coated  with  PDMAA-BP. 
The  intrusion  is  measured  in  situ  with  a  Tesa-Visio  300  video  measuring  instrument. 
The  minimum  intrusion  of  20  p,m  occurs  directly  along  the  channel  walls,  the  max¬ 
imum  intrusion  is  measured  in  the  middle  of  the  channel.  Error  bars  indicate  one 
standard  deviation.  The  error  is  large  due  to  surface  roughness  of  the  diffusion  layer. 


The  cross-sectional  area  of  the  uncoated  channel  is  determined 
to  be  0.95  mm2,  the  cross-sectional  area  of  the  channel  coated  with 
PDMAA-BP  amounts  to  1.02  mm2  (Fig.  3). 

Part  of  the  fuel  stream  passes  through  the  diffusion  layer  and 
thus  further  reduces  the  fluid  velocity  in  the  flow  field  channel.  This 
flow  rate  Qdl  can  be  estimated  by  rewriting  Darcy’s  Law,  where  r\ 
is  the  dynamic  viscosity  of  the  2-molar  methanol  solution,  l  is  the 
channel  length  touching  the  active  area  and  thus  the  diffusion  layer, 
and  Wchan  is  the  channel  width,  k  is  the  in-plane  permeability  [30] 
and  d  the  thickness  of  the  employed  Toray  TP-060  diffusion  layer. 
A p  is  the  pressure  loss  over  the  channel  length.  In  channels  with 
hydrophobic  walls,  this  is  mainly  dominated  by  the  capillary  pres¬ 
sure  loss  over  the  blocking  bubble,  so  we  assume  that  A p  =  A pbub: 

Qdl  =  — C^n-Apbub  (!) 


Thus,  the  average  fuel  velocity  vChan  in  the  employed  flow  field 
channel  can  be  calculated  as 


Vchan 


Qsp  -  Qdl 
Achan 


(2) 


with  the  flow  rate  generated  by  the  syringe  pump  Qsp  and  the  open 
channel  cross-sectional  area  ACftan.  In  the  channel  with  the  PDMA- 
PP  coating,  A p  is  close  to  zero  and  thus  we  expect  only  negligible 
diffusion  layer  flow. 

The  syringe  pump  is  calibrated  to  generate  a  fuel  flow  rate 
between  0  and  15.6 puls-1  which  leads  to  average  fuel  velocities 
of  0-1 6.3  mm  s-1  in  the  uncoated  channel  and  0-1 5.3  mm  s-1  mm 
in  the  coated  channel.  These  are  typical  velocities  for  DMFCs,  as  is 
verified  by  employing  Faraday’s  Law: 


dlCH3OH  = 


Mch3ohJ 

zF 


(3) 


where  the  mass  flow  of  methanol  mCH30H  consumed  by  the  reac¬ 
tion  is  calculated  from  the  generated  current  /,  the  molar  mass  of 
methanol  MCh3oh.  the  Faraday  constant  F  and  the  valency  number 
z.  For  example,  using  a  single-channel  meander  flow  field  with  a 
width  of  1  mm  and  a  height  of  1  mm,  a  typical  current  density  of 
100  mA  cm-2  [31]  and  a  2-molar  methanol  solution  pumped  at  a 
typical  stoichiometry  of  5  [32]  through  an  active  area  of  25  cm2 
would  lead  to  an  average  fuel  velocity  of  1 0.8  mm  s-1 . 

The  pressure  drop  over  a  single  blocking  bubble  is  calculated 
using  the  Young-Laplace  equation  for  the  difference  in  pressure  pc 
across  the  interface  between  two  different  fluids  with  the  interfa¬ 
cial  tension  a  and  the  local  curvature  K: 


Pc  =  orK  (4) 

The  pressure  drop  of  a  bubble  can  be  subsequently  formulated 
as: 


APbub  ~  Pc,f  Pc,b  (3) 

pCf  and  Pe  t,  describe  the  difference  in  pressure  across  the  liquid/gas 
interface  at  the  front  and  the  back  of  the  bubble,  respectively.  These 
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Fig.  4.  Schematic  diagram  of  a  pressure  drop  over  a  blocking  bubble,  which  is 
derived  by  subtraction  of  the  pressure  drop  across  the  meniscus  at  the  front  of  the 
bubble  from  the  pressure  drop  across  the  meniscus  at  the  back  of  the  bubble.  These 
are  in  turn  calculated  from  contact  angles  and  the  bubble  width  and  height. 


can  be  derived  from  the  principal  curvatures  between  the  left  chan¬ 
nel  wall  and  the  right  channel  wall  and  between  the  channel  floor 
and  the  diffusion  layer: 


Pcj  =  <r 

Pc,b  =  ^ 


f  COS(0rec,w)  +  COS(0rec  qi)  ^  2  COS ($rec,w) 
\  hub  Wbub 

f  cos(flad,w)  +  cos(flad;DL)  2cos(6>ad;W)> 

\  hub  wbub  ) 


(6) 

(7) 


with  the  height  hbub  and  width  wbub  of  the  bubble  in  the  channel, 
and  advancing  (ad)  and  receding  (rec)  contact  angles  0  on  the  wall 
(w)  and  the  diffusion  layer  (DL)  (Fig.  4). 

Contact  angles  of  the  flow  field  are  determined  using  the  sessile 
drop  method.  A  droplet  of  the  methanol  solution  used  in  the  DMFC 
was  applied  to  both  COC  with  a  PDMAA-BP  coating  and  native, 
uncoated  COC.  Advancing  (pre-wetted)  and  receding  contact  angles 
were  determined.  On  the  native  COC,  an  advancing  contact  angle  of 
68°  ±  3°  and  a  receding  contact  angle  of  45°  ±  3°  were  measured.  On 
the  COC  coated  with  PDMAA-BP,  the  advancing  contact  angle  was 
determined  to  be  10°  ±  3°,  the  receding  contact  angle  was  close  to 
0°.  In  addition,  both  advancing  and  receding  contact  angles  of  the 
diffusion  layer  are  estimated  to  be  close  to  0°,  as  the  experiments 
in  the  hydrophilic  flow  field  coated  with  PDMAA-BP  show  virtually 
no  pinning  of  blocking  bubbles. 

The  contact  angles  used  for  calculation  in  this  work  are  assumed 
to  be  independent  of  bubble  velocity.  At  higher  bubble  velocities, 
dynamic  contact  angles  would  start  to  play  a  role  and  initially, 
advancing  contact  angles  would  increase  [33].  The  bubble  shape 
would  change  as  the  velocity  rises  further,  which  would  lead  to  an 
increase  of  bubble  pressure  drop  and  in  consequence  to  an  increase 
of  relative  velocity  between  the  bubble  and  the  fuel  stream. 

We  conclude  that  dynamic  contact  angles  can  be  neglected  in 
the  reported  work.  This  is  because  we  observe  no  significant  shape 
change  once  the  bubble  is  in  motion  (Fig.  5).  We  estimate  the 
dynamic  contact  angle  for  a  bubble  velocity  of  1 5  mm  s-1  using  the 
formula  introduced  by  Bracke  et  al.  [34]: 

COS  ^ad,dyn  —  C°S  @ad,stat  —  2(1  +  COS  0ad,stat)^°'5  (8) 


with 


Ca  = 


PV bub 
(7 


(9) 


where  the  advancing  dynamic  contact  angle  0addyn  is  calculated 
from  the  corresponding  static  contact  angle  0adstat  and  the  dimen¬ 
sionless  velocity  expressed  as  the  capillary  number  Ca,  which  in 
turn  is  derived  from  the  dynamic  viscosity  p  and  surface  tension  a 
of  the  liquid  and  the  bubble  velocity  vbub.  For  the  maximum  bub¬ 
ble  velocity  observed  in  this  experiment,  this  would  increase  the 
advancing  contact  angle  on  the  channel  wall  Oad.w  by  3°  from  68° 
to  71°,  which  was  considered  to  be  negligible. 


3.  Results 


3.1.  Pinning  and  bubble  velocity  in  hydrophilic  PDMAA-BP  coated 
COC  flow  field 

Three  different  blocking  bubbles  of  a  total  length  lbub  of  1.4,  2.8 
and  4.8  mm  (Fig.  6)  were  generated  and  the  fuel  velocity  was  varied 
from  0  to  1 5.3  mm  s-1 .  We  observe  little  pinning  for  all  bubbles.  As 
soon  as  the  fuel  feed  is  activated,  the  bubbles  start  to  move.  As  we 
increase  the  flow  rate,  the  bubbles  equally  increase  their  velocities 
to  match  the  average  fuel  velocity. 

The  bubble  geometry  is  symmetric,  with  all  advancing  and 
receding  contact  angles  being  close  to  0°,  so  that  both  the  front 
of  the  bubble,  in  the  direction  of  flow,  and  the  back  of  the  bubble, 
facing  the  flow,  form  a  hemisphere  with  the  triple  line  between  the 
channel  floor,  methanol  solution  and  CO2  and  the  corresponding 
triple  line  between  the  diffusion  layer,  methanol  solution  and  CO2 
moving  in  parallel.  Small  fuel  bypass  channels  around  a  bubble  are 
formed  in  all  four  corners,  in  accordance  with  the  Concus  Finn  con¬ 
dition  [35].  This  states  that  the  sum  of  both  receding  contact  angles 
of  the  surfaces  forming  a  corner  and  the  corner  angle  is  not  allowed 
to  be  180°  or  greater  if  a  bypass  in  the  corner  of  a  channel  is  to  be 
formed.  Considering  the  two  different  surface  types  of  channel  and 
diffusion  layer,  we  can  formulate  the  two  respective  cases: 

@rec,DL  +  $rec,w  +  Oi  <  180°  (10) 

2  0rec,w  +  a<  180°  (11) 

In  this  experiment,  all  receding  angles  on  both  the  diffusion 
layer  and  the  channel  wall  and  floor  are  close  to  0°.  With  the  cor¬ 
ner  angle  of  90°  between  channel  wall  and  floor  and  an  angle  of 
approximately  83°  between  diffusion  layer  and  channel  wall,  due 
to  diffusion  layer  intrusion  into  the  channel,  all  sums  are  clearly 
lower  than  180°  and  bypass  formation  is  inevitable. 

On  plotting  bubble  velocity  versus  average  fuel  velocity  (Fig.  7), 
we  observe  a  slope  of  1 .02  for  all  bubble  sizes,  which  means  that  the 
bubbles  travel  slightly  faster  than  the  average  fuel  velocity.  This  can 
be  explained  by  the  parabolic  flow  profile.  Due  to  the  bypasses  to 
the  bubble  in  all  four  corners  of  the  channel,  where  the  fuel  flow  is 
close  to  the  walls  and  thus  flow  velocity  is  considerably  lower  than 
in  the  channel  centre,  the  bubble  is  in  contact  with  a  flow  region  of 
higher  average  velocity  than  the  total  average  channel  velocity. 


Fig.  5.  Bubble  geometry  at  two  different  average  flow  velocities  inside  a  native  hydrophobic  COC  channel:  (a)  2.6  mm  s-2,  (b)13.6  mm  s~2.  The  change  in  the  bubble  geometry 
is  negligibly  small  at  the  investigated  fluid  velocities.  The  bubble  is  depicted  schematically  in  Fig.  6g. 
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Fig.  6.  Bubble  observed  in  experiment  at  an  average  flow  velocity  of  5  mm  s-2 :  (a)-(c)  hydrophobic  native  COC  channel,  (d)-(f)  COC  channel  coated  with  hydrophilic  PDMAA- 
BP,  (g)  schematic  diagram  of  bubble  in  hydrophobic  native  COC  channel,  (h)  schematic  diagram  of  bubble  in  COC  channel  coated  with  hydrophilic  PDMAA-BP.  Bubble  length: 
(a)  2.4  mm  (b)  3.2  mm  (c)  5.4  mm  (d)  1.4  mm  (e)  2.8  mm  (f)  4.8  mm. 


Due  to  negligible  contact  angle  hysteresis,  we  do  not  observe 
any  change  in  the  shape  of  the  bubble  when  the  pumping  direction 
is  reversed. 

3.2.  Bubble  pinning  and  velocity  in  hydrophobic  native  COC  flow 
field 

Three  different  blocking  bubbles  with  a  total  length  lbub  of  2.4, 
3.2  and  5.4  mm  (Fig.  6)  were  generated  and  the  bubble  velocity 
was  plotted  versus  the  average  fluid  velocity  (Fig.  8).  Observing 
the  smallest,  2.4  mm  bubble,  we  identify  three  different  types  of 
behaviour  depending  on  fuel  velocity: 

At  a  low  velocity  of  up  to  2  mm  s-1  we  observe  pinning,  with  the 
contact  angle  hysteresis  increasing  at  the  walls  of  the  channel  as  the 
velocity  rises.  At  the  front  of  the  bubble,  in  the  direction  of  flow,  the 
receding  contact  angle  is  established.  At  the  back  of  the  bubble,  the 
advancing  contact  angle  is  established.  In  contrast  to  the  diffusion 
layer,  where  contact  angle  hysteresis  is  negligible,  the  hysteresis 
here  causes  the  contact  angle  at  the  channel  floor  to  differ  signifi¬ 
cantly  between  the  back  and  the  front  of  the  bubble.  Consequently, 
the  shape  of  the  bubble  caps  is  different  as  illustrated  in  Fig.  6.  At 
the  back  of  the  bubble  the  part  of  the  meniscus  touching  the  chan¬ 
nel  floor  is  pulled  back  away  from  the  bubble  centre  so  that  the 
receding  contact  angles  of  both  walls  and  the  diffusion  layer  can  be 


Fig.  7.  Bubble  velocity  versus  average  flow  velocity  for  3  bubble  lengths  in  the  COC 
channel  coated  with  hydrophilic  PDMAA-BP.  We  do  not  observe  any  pinning.  The 
bubble  velocity  is  independent  of  bubble  length. 


formed.  By  contrast,  the  difference  is  considerably  lower  at  the  front 
of  the  bubble  so  that  a  cap  is  formed,  which  more  closely  resembles 
a  hemisphere,  with  the  triple  line  between  the  channel  boundary, 
methanol  solution  and  CO2  on  the  channel  floor  slightly  further 
forward  than  the  triple  line  on  the  diffusion  layer.  Due  to  bubble 
pinning,  the  methanol  solution  must  bypass  the  bubble  in  the  cor¬ 
ners  of  the  channel.  We  do  not  observe  any  bypass  in  the  corners 
between  the  channel  walls  and  floor.  Instead,  the  fuel  overtaking 
the  bubble  bypasses  the  bubble  in  the  corners  between  the  channel 
walls  and  the  diffusion  layer,  leading  to  a  detachment  of  most  of  the 
meniscus  facing  the  diffusion  layer.  This  effect  is  further  reinforced 
due  to  the  fact  that  the  diffusion  layer  arches  into  the  channel,  as 
illustrated  in  Fig.  2,  so  that  the  bubble  only  rides  on  a  narrow  strip 
in  the  middle  of  the  diffusion  layer.  Therefore  the  main  part  of  the 
diffusion  layer  is  in  direct  contact  with  the  methanol  solution.  With 
the  contact  angles  measured  in  Section  2.2,  this  satisfies  the  Con- 
cus  Finn  condition  already  described  in  Eqs.  (10)  and  (11),  where 
the  sum  of  both  receding  contact  angles  of  the  surfaces  forming  a 
corner  and  the  corner  angle  is  not  allowed  to  be  180°  or  larger  so 
that  a  bubble  bypass  in  the  corner  of  a  channel  can  be  formed. 

At  2  mm  s-1 ,  the  bubble  has  a  fully  developed  contact  angle  hys¬ 
teresis  and  will  retain  this  shape  even  at  higher  velocities  as  long  as 
the  direction  of  the  flow  does  not  change.  At  a  higher  velocity  than 
2  mms-1,  the  bubble  detaches  itself  from  its  pinning  position  and 
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Fig.  8.  Bubble  velocity  versus  average  flow  velocity  for  3  bubble  lengths  in  the 

hydrophobic  native  COC  channel.  We  observe  pinning  in  dependence  on  the  bubble 
length.  Longer  bubbles  require  a  lower  average  flow  velocity  to  detach  themselves 
from  their  pinning  position. 
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starts  to  move  in  the  direction  of  the  flow.  Additionally,  we  observe 
that  between  2  and  3  mm  s-1 ,  the  bubble  is  subject  to  stop-and-go 
motion,  which  is  attributed  to  diffusion  layer  roughness.  The  bub¬ 
ble  movement  is  assumed  to  be  interrupted  by  carbon  fibres  that 
protrude  into  the  channel. 

Above  3  mm  s-1 ,  we  observe  a  linear  increase  of  bubble  velocity 
with  the  average  fluid  velocity.  The  gradient  is  slightly  less  than  one, 
so  the  relative  velocity  between  the  bubble  and  the  fluid  decreases 
slightly  in  the  velocity  range  investigated  in  this  work,  once  the 
bubble  is  in  motion. 

The  medium-sized  bubble  starts  moving  at  considerably  lower 
velocities.  We  observe  pinning  between  0  and  0.7  mm  s-1  fuel 
velocity.  There  is  no  stop-and-go  behaviour,  the  larger  bubble  copes 
better  with  the  surface  roughness  but  at  0.7  mm  s-1  and  faster  we 
observe  a  linear  rise  of  approximately  one. 

The  largest  bubble  generated  starts  to  move  at  even  lower 
velocities.  We  observe  pinning  merely  between  0  and  0.5mms_1 
fuel  velocity.  There  is  no  stop-and-go  behaviour.  Starting  from 
0.5  mm  s-1  we  observe  a  linear  rise  of  bubble  velocity  with  fluid 
velocity,  with  a  slope  slightly  greater  than  one. 

It  is  noteworthy  that  the  largest  bubble,  starting  at  about 
1 5  mm  s-1 ,  travels  faster  than  the  average  channel  velocity,  a  phe¬ 
nomenon  observed  and  explained  in  Section  3.1 .  Small  and  medium 
bubbles  do  not  show  this  behaviour. 

Additionally,  the  pumping  direction  was  frequently  reversed  to 
observe  the  change  in  bubble  geometry.  As  one  can  expect,  advanc¬ 
ing  and  receding  contact  angles  trade  locations  and  the  bubble 
moves  in  the  opposite  direction  with  the  same  absolute  value  of 
relative  velocity. 

3.3.  Impact  on  mass  transport  and  DMFC  performance 

The  pressure  gradient,  that  is  created  by  the  flow  to  develop 
and  maintain  the  full  contact  angle  hysteresis  against  friction  [36], 
and  the  ability  to  open  bypasses  along  the  bubble,  lead  to  a  relative 
velocity  between  the  gas  bubble  and  the  methanol  solution,  even 
at  fuel  velocities  that  allow  no  bubble  pinning.  The  decisive  factor 
determining  how  quickly  C02  can  be  transported  out  of  the  flow 
field  in  any  given  channel  and  diffusion  layer  configuration  is  the 
length  of  the  blocking  bubble.  This  configuration  can  in  turn  be  var¬ 
ied  by  using  different  materials  to  reduce  or  increase  the  contact 
angle.  Lower  contact  angles  will  not  only  facilitate  bubble  trans¬ 
portation  but  also  reduce  the  pressure  drop  per  blocking  bubble. 
Consequently,  the  energy  lost  by  the  pumping  needed  to  overcome 
the  flow  field  pressure  drop  can  be  decreased  by  reducing  the  con¬ 
tact  angles  of  the  boundaries  enclosing  the  channel.  At  first  sight, 
this  would  identify  the  channel  with  the  hydrophilic  PDMAA-BP 
coating  as  the  favoured  configuration  to  generate  higher  power 
densities. 

On  second  thought,  the  native  COC  channel,  featuring  a  highly 
hydrophilic  diffusion  layer  and  a  hydrophobic  channel,  offers  a 
unique  advantage  compared  to  the  completely  hydrophilic  con¬ 
figuration.  Due  to  the  fact  that  this  configuration  yields  bypasses 
only  in  the  corners  between  the  diffusion  layer  and  the  channel 
walls  in  combination  with  a  high  pressure  drop  between  the  back 
and  front  of  the  bubble,  high  fuel  velocities  can  be  expected  above 
the  diffusion  layer  under  a  blocking  bubble.  This  aids  methanol 
transport  to  and  reaction  product  removal  away  from  the  reaction 
centre,  so  that  mass  transfer  is  enhanced.  It  is  additionally  facil¬ 
itated  by  the  diffusion  layer  arching  into  the  channel  due  to  the 
clamping  pressure,  so  that  most  of  the  diffusion  layer  is  directly 
exposed  to  the  fuel  stream.  This  effect  increases  as  the  blocking 
bubble  density  increases,  which  offers  the  additional  advantage  of 
enhancing  mass  transport  in  dependence  on  channel  length,  which 
compensates  for  methanol  concentration  losses  due  to  reaction  and 
crossover. 


In  contrast  to  this,  bubble  transportation  in  completely 
hydrophilic  channels  will  only  cause  a  minor  disturbance  of  the  typ¬ 
ical  parabolic  flow  profile  in  the  middle  of  the  channel  and  the  fuel 
layer  directly  above  the  diffusion  layer  is  mostly  stagnant.  Thus,  the 
supply  of  methanol  and  waste  discharge  mainly  relies  on  diffusion. 
To  study  this  further,  the  generated  velocities  in  the  uncoated  chan¬ 
nel  configurations  are  calculated  using  a  model  which  is  described 
in  the  following  section. 

4.  Model 

4.1.  Description 

The  model  presented  is  a  three-dimensional,  two-phase  model 
of  a  blocking  bubble  in  the  flow  field  channel  of  a  DMFC  for  the 
purpose  of  investigating  fuel  transport  along  the  boundary  between 
the  diffusion  layer  and  the  channel.  A  detailed  description  of  the 
following  phenomena  is  provided  in  this  section: 

(1 )  Average  velocity  of  methanol  solution  bypassing  the  bubbles  in 
the  corners  between  the  diffusion  layer  and  channel  wall. 

(2)  Bubble  velocity  resulting  from  surface  properties,  bubble  length 
and  mean  flow  rate  of  the  fuel  in  the  channel. 

The  main  assumptions  in  this  model  are  as  follows: 

(1)  The  model  is  defined  to  be  isothermal. 

(2)  All  fluids  are  incompressible. 

(3)  All  flow  is  laminar. 

(4)  Concentration  changes  due  to  chemical  reactions  are  neglected. 

(5)  The  bubble  is  considered  to  be  in  motion,  which  means  that 
the  contact  angle  hysteresis  is  fully  developed,  and  the  bubble 
shape  is  assumed  to  be  constant  and  independent  of  transport 
velocity. 

(6)  The  channel  floor  and  walls  are  defined  to  be  no-slip  walls. 

(7)  Bubble  menisci  are  treated  as  slip  boundaries. 

(8)  Only  advancing  and  receding  contact  angles  are  used;  a  gradual 
transition  over  the  length  of  a  bypass  is  not  considered.  Conse¬ 
quently,  a  fully  developed  bypass  is  defined  as  a  channel  with 
a  constant  cross-section. 

4.2.  Model  domain,  equations  and  numerical  solution 

The  computational  domain  employed  for  the  simulation  part  of 
this  work  is  depicted  in  Fig.  9a.  The  domain  represents  a  section 
containing  one  blocking  bubble,  with  the  channel  cut  perpendicu¬ 
larly  to  the  diffusion  layer,  channel  walls  and  channel  floor  along 
the  front  and  back  of  the  bubble.  The  frame  of  reference  is  the  gas 
bubble  with  the  coordinate  system  fixed  at  the  centre  of  gravity 
of  the  bubble  so  that  the  simulation  of  bubble  motion  is  described 
by  moving  walls  and  corresponding  methanol  solution  flow.  The 
computational  domain  is  divided  into  subdomains  £2[  with  £2\ 
describing  the  parts  of  the  channel  around  the  bubble  filled  with 
methanol  solution  and  Q2  defined  as  the  gas  bubble.  These  sub- 
domains  are  coupled  with  each  other  and  enclosed  by  boundary 
conditions  The  inlet  and  outlet  of  the  domain  are  described 
by  <5^1  and  8Q2  respectively.  Both  are  defined  as  constant  pres¬ 
sure  boundary  conditions,  with  8£2\  set  to  the  sum  of  ambient 
pressure  and  the  pressure  drop  from  inlet  to  outlet.  8£22  is  set  to 
ambient  pressure.  8£23  describes  the  bubble  meniscus  between  the 
methanol  solution  and  C02  gas  and  is  defined  as  a  slip  boundary. 
At  8£24  to  8£27<  no-slip  boundary  conditions  between  fuel  and  wall 
[8Q4\  fuel  and  diffusion  layer  (&f25),  bubble  and  wall  (5-f26)  and 
bubble  and  diffusion  layer  (8£27 )  are  set. 
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Fig.  9.  (a)  Model  domain:  liquid  phase  Q2 :  gas  phase;  interfaces:  8Q\\ 

methanol  solution  inlet  8£22'-  methanol  solution  outlet  8P2y.  bubble  meniscus 
between  methanol  solution  and  CO2  gas  8 £2 4:  methanol  solution channel  wall 
8P2 5:  methanol  solutions  diffusion  layer  8P26\  C02  gas  channel  wall  8P2 7:  CO2 
gas  diffusion  layer  and  (b)  we  assume  that  the  capillary  pressure  drop  over 
the  gas  bubble  equals  the  pressure  drop  over  the  bubble  bypasses  and  can  thus 
calculate  the  volume  flow  by  reconstructing  the  bypass  and  using  it  as  our  simpli¬ 
fied  model  domain.  Model  domain:  Q\\  liquid  phase;  interfaces:  8P2\\  methanol 
solution  inlet  8£22\  methanol  solution  outlet,  (c)  Cross-section  of  the  domain 
described  in  (b)).  Model  domain:  Q\ :  liquid  phase;  interfaces:  8P2 3:  bubble  menis¬ 
cus  between  methanol  solution  and  C02  gas  8P2 4:  methanol  solution  channel  wall 
8P2  5:  methanol  solution  ^diffusion  layer;  contact  angles:  0W:  methanol  solution  on 
COC  wall  9dl  :  methanol  solution  on  diffusion  layer. 


To  reduce  calculation  time  and  allow  for  parametric  studies,  the 
model  is  further  simplified  by  the  following  assumptions: 

The  pressure  distribution  in  the  gas  within  the  bubble  can  be 
considered  to  be  constant. 

The  capillary  pressure  drop  over  the  bubble,  which  has  to  be 
matched  by  the  liquid  phase  of  the  channel,  will  mainly  occur 
in  the  bypasses  along  the  bubble  due  to  their  relatively  small 
cross-section  compared  to  the  cross-section  of  the  channel.  So  we 
assume: 

APbub  —  APbyp  (12) 

with  the  capillary  pressure  drop  over  the  bubble  A pbub  and  the 
linear  pressure  drop  A pbyp  in  the  bypass  channel  of  the  bubble. 

The  slip  boundary  can  be  approximated  by  a  symmetry  condi¬ 
tion  which  does  not  allow  any  momentum  transfer. 

These  assumptions  lead  to  a  simplified  domain  which  is  depicted 
in  Fig.  9b  and  c,  where  a  single  bypass  containing  methanol  solu¬ 
tion  features  constant-pressure  boundaries  at  the  inlet  (8£2\ )  and 
outlet  (8  £2 2),  and  symmetric  conditions  (<5J23)  describe  the  bubble 
meniscus,  as  well  as  the  no-slip  walls  of  the  channel  (8 £2 4)  and  the 
diffusion  layer  {8  £2  5). 

All  densities  are  considered  constant  so  mass  conservation  may 
be  reduced  to: 

VpF  =  0  (13) 

with  the  velocity  vector  v  and  the  density  of  methanol  solution  p. 

In  general,  the  velocity  field  is  determined  by  the  Navier-Stokes 
equation  for  incompressible  fluids: 


non-zero.  So  we  conclude:  v  •  Vv  =  0.  In  analogy  to  this,  the  pres¬ 
sure  similarly  depends  only  on  z  and  is  thus  a  linear  function  of 
z.  Due  to  the  fact  that  the  coordinate  system  is  fixed  at  the  cen¬ 
tre  of  gravity  of  the  bubble  and  we  hence  only  study  a  progression 
of  velocity  fields  in  steady  state,  we  can  additionally  assume  that 
dv/dt  equals  zero.  We  can  now  formulate  what  is  commonly  known 
as  Hagen-Poiseuille  flow  [37]: 

(9*  +  dy)vz(x,y)  =  (15) 

with  the  length  l  of  the  channel  and  the  pressure  drop  A  p  from  inlet 
to  outlet. 

This  expression  for  the  velocity  field  can  now  be  used  to  derive 
the  volume  flow  rate  in  the  channel  as  a  function  of  the  pressure 
drop  for  an  arbitrary  cross-sectional  shape  by  the  introducing  a 
specific  correction  coefficient  C  [37]: 


 ApA2 


(16) 


As  both  the  volume  flow  lv  and  the  correction  coefficient  C  are 
unknown,  we  apply  CFD  to  solve  this  problem  numerically,  as 
described  in  the  following  section. 


4.3.  Model  validation:  comparison  of  modelling  results 

The  geometrical  configuration  of  the  bubble  bypass  channels  is 
reconstructed  from  the  films  of  the  experiments.  The  bypass  chan¬ 
nels  are  described  by  the  bypass  length  lcor  and  bypass  height  hcor. 
With  the  known  receding  contact  angles  at  the  wall  6rec,w  and  dif¬ 
fusion  layer  0rec,Dh  the  cross-sectional  areaAcor  is  well-defined.  The 
bypass  is  modelled  in  the  CAD  SolidWorks  software,  imported  into 
the  CFD-GEOM  meshing  tool  and  meshed  with  a  structured  grid 
consisting  of  50,000  to  100,000  cells  depending  on  the  size  of  the 
bubble  bypass  channels.  We  assume  that  A p  is  equal  to  the  cap¬ 
illary  pressure  drop  Apbub  over  the  bubble,  which  in  turn  can  be 
calculated  as  demonstrated  in  Section  2.2. 

Finally  the  model  is  solved  by  applying  the  CFD-ACE  +  simulation 
software  from  the  ESI  Group  company  [38].  As  the  bubble  serves  as 
the  frame  of  reference,  the  bubble  velocity  is  described  indirectly  by 
moving  the  walls  (8 £2 4, 8  £2$)  against  the  direction  of  flow  with  the 
desired  magnitude  of  velocity.  This  procedure  can  be  further  simpli¬ 
fied  by  only  simulating  the  transition  point  between  bubble  pinning 
and  movement,  where  the  bubble  velocity  vbub  is  still  zero  but  the 
contact  angle  hysteresis  is  fully  developed.  After  convergence,  a 
bubble  bypass  flow  rate  Iv  is  calculated  and  from  this,  a  correc¬ 
tion  coefficient  C  for  the  specific  geometric  form  of  the  respective 
cross-section  is  derived: 

q  _  ^ Pbub^cor  (17) 

iylcorP 

The  starting  point  to  derive  the  relation  between  the  average 
fuel  velocity  in  the  channel  and  the  bubble  velocity  is  the  balance  of 
flow  rates  between  the  channel  (subscript  chan),  bypass  (subscript 
cor)  and  bubble  (subscript  bub)  at  the  transition  point  that  can  be 
described  as: 

^chan^chan  =  V bubi^chan  —  2ACor  )  +  2  VcorAcor  (18) 


p  +  v  •  VDJ  =  -Vp  +  pV2i>  (14) 

with  the  velocity  vector  v,  the  static  pressure  p  and  density  p  and 
viscosity  p  of  the  methanol  solution.  As  stated  above,  our  model 
domain  is  assumed  to  be  a  channel  parallel  to  the  z  axis  with  a 
constant  cross-section  in  the  xy-plane  (Fig.  9).  Consequently,  the 
forces  in  the  xy-plane  cancel  each  other,  leading  to  a  velocity  field 
independent  of  z  where  only  the  sum  of  its  z  components  can  be 


With 


K  —  VcorA 


cor 


(19) 


and  Eq.  (17),  vCOr  is  substituted  in  Eq.  (18)  and  a  general  dependence 
of  the  bubble  velocity  vbub  on  the  average  channel  velocity  vcban  can 
be  formulated: 


vbub  ■ 


^chan^chan 

\han  ~  2 ACoi 


2  A  pACOr^ 

Cplcor{Acban  —  2  ACOr) 


(20) 
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Fig.  10.  Comparison  of  simulation  and  experiment:  Bubble  velocity  versus  average 
channel  velocity  for  3  bubble  lengths  in  the  COC  channel  coated  with  hydrophilic 
PDMAA-BP.  Excellent  agreement  between  simulation  and  experiment  is  achieved. 


Fig.  12.  Simulation  results:  bubble  velocity  and  average  bubble  bypass  velocity 
versus  average  channel  velocity  for  3  bubble  lengths  in  the  hydrophobic  native  COC 
channel.  Shorter  bubbles  induce  considerably  higher  fuel  velocities  in  their  bypass 
channels. 


This  approach  is  repeated  to  cover  the  range  of  bubble  velocities 
and  geometric  dimensions  observed  in  the  experimental  part  of  this 
work.  Finally  the  results  are  compared  to  the  experimental  data  of 
both  coated  (Fig.  10)  and  uncoated  channel  (Fig.  11). 

Excellent  agreement  between  simulation  and  experiment  is 
achieved. 


4.4.  Average  fuel  velocity  bypassing  the  bubble 


We  reconsider  the  hypothesis  formulated  in  Section  3.3,  where 
high  velocities  of  methanol  solution  in  the  bubble  bypasses  of  the 
hydrophobic  channel  directly  above  the  diffusion  layer  were  pre¬ 
dicted.  We  can  now  derive  a  formula  to  predict  the  fuel  velocity 
in  the  bypass  of  a  blocking  bubble  relative  to  the  channel  wall  by 
replacing  vbub  in  Eq.  (18)  with  Eq.  (20): 


_  A pAcor 
Cplcor 


(21) 


In  conclusion,  this  means  that  the  relative  velocity  of  the  fuel 
bypassing  a  blocking  but  non-pinning  bubble  is  independent  of 
both  bubble  and  mean  fuel  velocity  inside  the  channel,  as  soon  as 
contact  angle  hysteresis  and  subsequently  the  pressure  loss  over 
the  bubble  is  fully  developed.  Blocking  bubbles  in  hydrophilic  chan¬ 
nels  with  a  pressure  loss  of  approximately  zero  will  consequently 
cause  a  momentary  standstill  in  fuel  transport.  Increasing  the  con¬ 
tact  angle  hysteresis  will  lead  to  higher  fuel  velocities  in  this  critical 
zone  directly  above  the  diffusion  layer. 

By  analysing  the  experimentally  generated  bubbles,  we  can  on 
the  one  hand  conclude  that  both  viscosity  p,  due  to  negligible 
changes  in  temperature  and  substance  concentration,  and  the  geo- 
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Fig.  11.  Comparison  of  simulation  and  experiment:  Bubble  velocity  versus  aver¬ 
age  channel  velocity  for  3  bubble  lengths  in  the  hydrophobic  native  COC  channel. 
Excellent  agreement  between  simulation  and  experiment  is  achieved. 


metric  correction  coefficient  C,  due  to  lack  of  significant  change 
in  cross-section  form,  can  be  regarded  as  constant  values  for  any 
given  bubble  length  in  our  system.  On  the  other  hand,  the  ratio  of 
the  bypass  cross-sectional  area  Acor  to  the  bypass  length  /cor  signif¬ 
icantly  decreases  as  the  blocking  bubble  increases  in  length,  and 
subsequently  the  fuel  velocity  vCOr  bypassing  the  bubble  decreases 
by  the  same  ratio.  In  the  native  hydrophobic  COC  channel  analysed 
in  this  work,  the  2.4  mm  long  bubble  induces  a  very  high  average 
velocity  of  74  mm  s-1  of  the  fuel  in  the  bypass  channel.  The  velocity 
decreases  as  the  blocking  bubble  becomes  longer.  The  3.2  mm  bub¬ 
ble  still  offers  an  average  velocity  of  the  fuel  of  30  mm  s-1  when  the 
bubble  is  bypassed.  Finally,  at  5.4  mm  bubble  length,  a  fuel  velocity 
of  9  mm  s-1  can  be  achieved  (Fig.  12). 

5.  Conclusion 

We  visually  investigated  blocking  CO2  bubbles  as  they  occur 
in  the  anode  of  a  direct  methanol  fuel  cell  (DMFC).  For  the  first 
time,  we  provide  a  side  view  of  the  two-phase  flow  in  the  chan¬ 
nel  by  fabricating  an  optically  transparent  flow  field  made  of  cyclic 
olefin  copolymer  (COC).  We  introduce  an  in  situ  method  to  quantify 
clamping-induced  flow  channel  intrusion  of  the  diffusion  layer  by 
employing  a  Tesa-Visio  300  video  measuring  instrument.  A  maxi¬ 
mum  intrusion  of  50  p,m  in  the  middle  of  the  channel  is  measured. 

We  compare  blocking  bubble  behaviour  in  a  native,  hydropho¬ 
bic  COC  channel  to  a  channel  coated  with  hydrophilic  PDMAA-BP 
at  average  methanol  solution  velocities  of  up  to  16.3  mms-1.  We 
use  the  same  hydrophilic  diffusion  layer  in  both  cases.  The  bub¬ 
ble  in  the  hydrophilic  channel  shows  almost  no  pinning  and  thus 
travels  at  approximately  the  same  mean  velocity  as  the  fuel.  In  con¬ 
trast  to  this,  a  blocking  bubble  in  the  hydrophobic  channel  pins  to 
the  channel  wall  at  low  mean  fuel  velocities  and  requires  consid¬ 
erable  fuel  velocities  to  reach  its  tear-off  point,  which  depends  on 
the  bubble  length.  The  critical  fuel  velocities  were  in  the  range  of 
0.5  mm  s-1  for  the  longest  and  2  mm  s-1  for  the  shortest  bubble 
studied  in  this  work.  The  receding  contact  angle  at  the  front  and  the 
advancing  contact  angle  at  the  back  of  a  blocking  bubble  induce  a 
significant  pressure  gradient  along  the  bubble.  This  gradient  leads 
to  a  relative  velocity  between  mean  fuel  velocity  and  blocking  bub¬ 
ble  which  is  approximately  constant,  even  at  fuel  velocities  inside 
a  flow  channel  that  do  not  allow  any  bubble  pinning. 

To  improve  the  understanding  of  bubble  dynamics,  a 
CFD/analytical  model  is  developed  that  describes  fuel  trans¬ 
port  around  C02  bubbles.  The  model  is  validated  and  excellent 
agreement  with  the  experimentally  measured  bubble  velocity  is 
achieved.  We  use  the  model  to  quantify  the  average  fuel  velocity 
of  methanol  solution  bypassing  a  blocking  bubble.  This  amounts 
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to  74mms-1  for  the  smallest  bubble  in  the  hydrophobic  channel 
investigated  in  this  work,  which  is  considerably  higher  than  the 
average  fuel  velocity  in  the  flow  channel.  This  leads  to  the  following 
conclusions  when  designing  flow  fields  and  diffusion  layers  for 
actively  pumped  DMFCs. 

The  losses  in  these  systems  must  be  evaluated  carefully.  On  the 
one  hand,  a  significant  fraction  of  the  generated  energy  is  wasted 
on  operating  a  fuel  supply  pump  which  is  strong  enough  to  over¬ 
come  the  pressure  losses  of  the  two-phase  flow  in  the  fuel  supply 
channels.  In  some  cases  it  is  advisable  to  employ  a  hydrophilic  flow 
field  to  minimise  the  total  pressure  drop  und  thus  open  up  the 
possibility  to  downgrade  to  a  pump  with  lower  power  consump¬ 
tion.  On  the  other  hand,  it  can  be  more  favourable  to  employ  the 
bubble  bypass  configuration  introduced  in  this  work.  We  combine 
a  flow  channel  surface  with  a  receding  contact  angle  of  >45°  and 
a  diffusion  layer  with  a  receding  contact  angle  close  to  0°.  In  this 
configuration,  fuel  bypass  channels  are  automatically  generated  in 
a  favourable  position  underneath  a  blocking  bubble  next  to  the  dif¬ 
fusion  layer.  At  the  same  time,  no  bypass  channels  exist  above  the 
blocking  bubble  away  from  the  diffusion  layer.  This  forces  fuel  to 
bypass  the  bubble  at  a  high  relative  velocity  close  to  the  diffusion 
layer,  which  is  very  efficient  for  supplying  fresh  fuel  to  the  DMFC. 
The  pressure  drop,  which  drives  fuel  through  these  bypass  chan¬ 
nels  at  considerably  higher  speed  than  the  average  channel  velocity, 
results  from  the  contact  angle  hysteresis  at  the  channel  walls  and 
floor  and  significantly  enhances  mass  transfer  in  the  DMFC. 

In  recent  years,  passive  DMFC  systems  have  become  increas¬ 
ingly  popular.  Obviously,  these  systems  do  not  have  any  pumping 
losses,  but  they  are  limited  instead  by  a  mediocre  power  density 
due  to  high  mass  transfer  losses.  This  case  considerably  favours 
hydrophobic  walls  and  a  hydrophilic  diffusion  layer  to  direct  back- 
flow  along  the  diffusion  layer,  as  the  bubble  rises  due  to  buoyancy, 
thus  enhancing  mass  transfer. 
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